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PREDICTION  OF  DRIFT  FORCES  ON  TWIN  HULL 
BODIES  IN  WAVES 


Choung  M.  Lee 
and 

Yoon-Ho  Kim 


David  W.  Taylor  Naval  Ship  UNITED  STATES 

Research  and  Development  Center 
Bethesda,  Maryland  20084 


A  theoretical  method  is  presented  to  predict  the  steady  mean  surge 
and  sway  forces  and  roll  and  yaw  moments  on  floating  twin-hull  ships 
in  waves.  The  steady  mean  forces  and  moments  are  derived  by  utiliz¬ 
ing  the  conservation  of  momentum  of  the  fluid  surrounding  the  body. 
The  associated  velocity  potentials  are  obtained  by  the  strip  approxi¬ 
mation.  Numerical  results  are  presented  in  graphs  and  pertinent 
discussions  are  made.  It  is  found  that  the  maximum  magnitude  of  the 
drift  force  on  SWATH  ships  would  be  greater  than  that  of  equivalent- 
displacement  monohull  ships  and  the  mean  heel  angle  induced  by  waves 
could  be  significant  for  SWATH  ships. 

1 .  Introduction 

It  is  a  well-known  physical  phenomenon  that  when  a  floating 
body  is  subject  to  surface  waves,  the  body  not  only  undergoes  an 
oscillatory  motion  but  also  its  mean  position  shifts  due  to  the 
drift  forces  and  moments  acting  on  the  body. 

The  drift  forces  are  often  referred  to  as  ’’second-order  mean 
forces.”  The  second-order  mean  forces  can  be  normally  determined  by 
utilizing  the  first-order  velocity  potential  solution  for  restrained 
and  unrestrained  bodies.  There  have  been  numerous  investigations 
conducted  on  the  second-order  mean  forces  such  as  added  resistance 
in  waves,  drift  force  and  moment,  and  suction  force  and  moment  for 
submerged  bodies.  These  investigations  have  been  applied  to  mono¬ 
hull  ships,  and  a  few  familiar  ones  are  given  in  References  [1-6). 

In  this  paper  an  analysis  of  the  second-order  mean  forces  and 
moments  in  the  horizontal  plane  is  carried  out  for  twin-hull  ships 
at  zero  forward  speed.  The  analysis  for  the  mean  surge  and  sway 
forces  and  yaw  moment  is  an  extension  of  the  work  of  Newman  [2] , 


and  the  analysis  for  the  steady  roll  moment  is  a  new  derivation. 
Both  analyses  are  based  on  the  principle  of  momentum  conservation 
of  the  fluid  surrounding  a  body. 

For  monohull  ships  the  steady  mean  heel  moment  induced  by 
waves  is  considered  to  be  insignificant  compared  to  the  steady  yaw 
moment,  since  there  exists  a  restoring  mechanism  for  the  roll  while 
there  is  no  such  mechanism  for  the  yaw.  The  foregoing  reasons 
could  have  been  the  cause  for  the  absence  of  the  derivation  for  the 
steady  mean  roll  moment.  However,  for  twin-hulls  such  as  a  semi- 
submersible  configuration  which  has  a  relatively  large  draft  as 
well  as  a  high  vertical  center  of  gravity,  it  is  not  obvious  if  the 
steady  mean  roll  moment  induced  by  waves  could  be  insignificant. 

One  of  the  main  purposes  of  this  paper  is  to  compute  the  steady 
mean  roll  moment  for  a  SWATH  ship  and  find  out  if  this  quantity 
should  be  treated  as  important  as  the  drift  force  and  the  steady 
mean  yaw  moment. 

The  solution  for  the  velocity  potential  function  involved  is 
obtained  by  a  strip  approximation  which  is  commonly  exercised  in 
the  prediction  of  ship  motion  in  waves  for  monohulls  [7]  and  twin- 
hulls  [81. 

Numerical  results  of  a  small-waterplane-area,  twin-hull 
(SWATH)  ship  are  presented  in  figures.  The  steady  mean  sway  force 
of  a  monohull  ship  is  also  computed,  and  the  results  are  compared 
with  those  of  the  SWATH  ship. 


2 .  Analysis  for  Mean  Forces  and  Moments 

The  fluid  is  assumed  incompressible,  inviscid,  and  its  motion 
irrota* ional .  The  depth  of  the  fluid  is  infinite,  and  a  chain  of 
progressive  plane  waves  represents  the  free  surface.  All  motions 
of  body  and  fluid  are  assumed  small  such  that  boundary  conditions 
can  be  linearized  and  imposed  on  the  undisturbed  positions  of  the 
body  and  on  the  calmwater  plane.  The  body  in  consideration  con¬ 
sists  of  twin-hulls  of  slender  configuration  which  are  rigidly 
connected  above  the  water.  The  slenderness  of  each  hull  is  such 
that  the  longitudinal  component  of  the  surface  normal  is  small 
compared  to  the  transverse  components.  The  body  has  no  forward 
speed  but  is  free  to  respond  to  the  incident  waves. 

Cartesian  coordinates  (x,y,z)  are  taken  such  that  the  origin 
is  located  on  the  calmwater  plane  z  *  0  directly  above  or  below  the 
center  of  gravity  of  the  body  at  rest,  the  positive  x-axis  is  di¬ 
rected  toward  the  bow,  and  the  positive  z-axis  is  vertically  upward 


(see  Figure  1).  The  fluid  motion  is  assumed  to  be  harmonic  in 
time  with  the  incident  wave  frequency. 


Figure  1  —  Description  of  Coordinate  System  and  Wave-Heading  Angle  (jj) 


The  fluid  velocity  V  at  any  point  can  be  represented  by  the 
gradient  of  a  velocity  potential  function  ♦  in  the  form 

V  =  V4>(x,y,z,t)  =  Re  [V  <J>  ( x ,  y ,  z )  e-  1Uit  ]  (1) 


where  Re  means  the  real  part  of  what  follows,  and  the  complex  func¬ 
tion  $  is  composed  of  the  incident  wave  potential  4>j.  and  the  dis¬ 
turbance  velocity  potential  <$>B,  i.e., 

$  ~  4>I  + 


The  incident  wave  potential  is  given  by 
igC» 

=  -  — —  exp(Kz  +  iKx  cos  \i  -  iKysiny)  (2) 

a 

where  g  is  the  gravitational  acceleration,  K  *  — ,  c*  the  wave 
amplitude,  and  y  the  angle  of  wave  incidence  relative  to  the  posi¬ 
tive  x-axis. 


The  disturbance  potential  $B  is  determined  by  the  solution  of 
the  Laplace  equation  satisfying  the  following  boundary  conditions: 
*Bz  '  k$b  *  0  *  on  z  -  0;  (3a) 


yBn 


In 


+  V 


(3b) 


*When  the  spatial  variables  x,  y,  z  and  the  temporal  variable  t 
are  used  as  a  subscript  of  a  function,  it  means  the  partial 
derivative  of  the  function  with  the  respective  variable. 
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on  the  body  surface  S0  at  its  rest  condition,  where  the  subscript  n 
indicates  a  normal  derivative  (the  unit  normal  vector  n  is  positive 
pointing  out  of  the  fluid  boundary),  and  VR  the  normal  velocity  of 
the  surface; 

■  0  as  z-»  -  •  (3c) 


Urn  /lT^  -  iK*B).  0  (3 

where  R  =  /(x-£)2  +  (y-n)2  anc*  (C,n)  represents  a  point  on  S0 . 

The  far-field  behavior  of  is  well  known  and  is  given  by 
Wehausen  and  Laitone  [8]  as 

/  v  /  w  v 


♦B  =  H(K,ir+0)exp|lCz  +  iKR  +  i 

where  H(k,8),  referred  to  as  the  Kochin  function,  is  defined  by 

H(k,0)  =JJ  exp(kc  +  xkC  cos  0  +  ikn  sin  0)ds  (4) 

S0  1 

From  the  conservation  principle  of  the  rate  of  change  of  the 
linear  and  angular  momentum  in  the  fluid  volume  V  enclosed  by  the 
surfaces  SQ ,  SF  (z  *  0  plane  excluding  the  waterplane  area  Sw) 
and  Sw  (the  far-field  control  surface),  we  can  obtain  the  expres¬ 
sions  for  the  mean  surge  force  F  ,  sway  force  F  and  yaw  moment 
—  *  1  2  — 

F6  (see  e.g.,  Newman  [2]  )  and  the  mean  roll  moment  F%  (see 


Appendix  I )  as 

v2  r  2* 

=  _  iyL  I  |H(K,ir+0)|*  cos  0  d0 

IT 

+ 

o»cA  cos  u  Im  H(K,  n-p) 

(5) 

«, 2  r  2n 

£{L  /  1  H(K ,  n+6  )  |  2  sin  0  d0 

**  J  0 

- 

^  u>CA  sin  y  Re  H(K,fr-u); 

(6) 

“  -S 

J-  Im  f**  H*(K,0)  ^  H(K,e)d6 

+ 

5a  Im  h  H(K’1T‘l,) 

(7) 

*The  slight  difference  in  the  expressions  are  due  to  the  differ¬ 
ence  in  the  expression  of  the  incident-wave  potential. 
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(8) 


where  the  asterisk  means  the  complex  conjugate  of  the  function, 

f  *  f  <IB>  +  f  <BB> 


*  H 

where 


¥ 


(IB)  ,  p<aCA 


Re^y  (zsiny  ♦  iy)  (♦ 
S* 


(^Bn  ^B  3n) 


aK(z  -  ix  cos  u  +  iy  sin  u)dg 


and 


^^(BB)  =  _  P*1  R eJ2U  d0  H*(Kfe)jy  (zsin6  +  iy) 


(*Bn  ~  *B*n) 


3  \  K(z  +  ix  cos  0  +  iy  sin  0). 
pBn  ‘  *B  5n  P 


e 

2r 


- 2-j  Re  i  f  d0^ 

16ir2  J  J 


dk 


k(k+K) 
"k-K  “ 


»*// 


(z  sin  0  +  iy ) 


(*Bn  ‘  *B  5s)  6 


0  0  5>0 

9_\  k(z  +  ix  cos  0  +  iysin6)ds 


(8a) 


(8b) 


It  should  be  noted  that  the  expressions  for  the  F^'s  shown  above  con¬ 
sist  of  two  terms  which  are  0(H2)  and  0(H).  From  Equations  (8a)  and 
(8b),  F„(IB)  =  0(H)  and  F„(BB)  =  0(H2).  We  can  readily  identify  that 
the  terms  of  0(H)  in  the  above  equations  represent  the  interaction  be¬ 
tween  ♦  j  and  4>fi,  and  the  terms  of  Q(H2)  represent  the  interaction  of 
<frB  by  itself. 

Within  the  linear  analysis  the  disturbance  potential  $B  can 
be  expressed  as  the  sum  of  the  diffracted  wave  potential  $7  and 

the  motion-related  potentials  for  J  -  1,  2 . 6,  i.e., 

7  J 

*B  ■  E  (•> 

where  Cj  for  j  -  1,  2 .  6  is  the  complex  amplitude  of  the 

j  th  mode  of  motion,  i.e.,  j  »  1,  2  and  3  correspond  to  surge,  sway 
and  heave,  respectively  and  J  ■  4,  5  and  6  correspond  to  roll,  pitch 
and  yaw,  respectively,  and  is  equal  to  the  incident  wave  ampli¬ 
tude 

For  the  solution  of  an  assumption  of  two-dimensional  flow 
condition  at  each  cross  section  is  made,  and  application  of  the 
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method  of  distributing  pulsating  sources,  which  satisfy  the  free- 
surface  condition  (3a),  on  the  immersed  contours  of  the  cross 
sections  is  made  [9,10].  In  this  way,  the  hydrodynamic  inter¬ 
actions  between  the  two  cross  sections  in  the  same  y-z  planes  are 
accounted  for  in  the  solution  of  .  The  motion  amplitudes  are 
obtained  from  coupled  linear  equations  of  motion  using  strip  theory 

[HI. 

The  Kochin  function  H  can  be  expressed  similarly  by 
7 

H  -  £  CjHj  (10) 


where 


Vk’0)  =//(* jn  "  *j  ?h) 


ek(c  +  i?  cos  0  +  insin0)dg 


3.  Results  and  Discussions 

A  twin-hull  ship  of  semi-submersible  configuration,  often  re¬ 
ferred  to  as  small-waterplane  area,  twin-hull  (SWATH),  is  chosen  for 
a  sample  calculation.  Its  particular  dimensions  are  presented  in 
Table  1  together  with  a  monohull  ship.  The  drift  force  computed  for 
this  monohull  will  be  compared  later  with  that  of  SWATH. 

As  can  be  seen  from  the  expressions  for  the  steady  mean  forces 
and  moments  given  by  Equations  (5)  through  (8),  the  Kochin  function 
H  is  directly  involved  in  these  formulae.  To  obtain  H  the  disturb¬ 
ance  potential  should  be  known  as  shown  by  Equation  (4).  The 
potential  is  further  decomposed  into  seven  components  which  are 
associated  with  the  six  modes  of  motion  and  the  wave  diffraction  as 
indicated  by  Equation  (9).  Thus,  the  necessary  computations  should 
begin  with  obtaining  the  for  j  *  1,  2,  .  .  .  ,7.  The  solutions 

of  ♦j's  are  utilized  to  obtain  the  hydrodynamic  coefficients  such  as 
the  added  mass  and  damping  coefficients  and  wave  excited  forces  and 
moments.  It  is  then  followed  by  the  computation  of  motion  in  six- 
degrees  of  freedom.  The  details  leading  to  this  stage  are  described 
in  Lee  and  Curphey  [11] ,  and  will  not  be  repeated  here. 


The  extra  condition  imposed  on  i.e.,  -  0  on  S g  (see  Appendix 

I)  is  carried  out  in  the  numerical  evaluation  of  in  this  work. 
This  condition  is  also  used  to  remove  the  so-called  irregular  fre¬ 
quencies  pointed  out  by  John  [12]. 
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Table  1 

Particular  Dimensions  of  SWATH  and  Monohull  in  Full  Scale 


Parameter 

SWATH 

Monohull 

Displacement,  long  tons  in  s.w. 

2802.0 

1761.0 

Length  at  the  Waterline,  m 

S2.S 

- 

Length  of  Main  Hull,  m 

73.2 

73.2 

Beam  of  Each  Hull  at  the  Waterline, m 

2.2 

13.4 

Hull  Spacing  Between  the  Centerlines,  m 

22.9 

- 

Draft  at  the  Midship,  m 

8.1 

3.62 

Bridging  Structure  Clearance  from 
Waterline,  m 

6.1 

_ 

Maximum  Diameter  of  Main  Hull,  m 

4.6 

- 

Longitudinal  Center  of  Gravity  Aft  of  Main 
Hull  Nose,  m 

36.5 

38.0 

Vertical  Center  of  Gravity  above 

Baseline,  m 

10.4 

5.3 

Transverse  GW,  m 

2.9 

1.91 

Longitudinal  5M,  m 

6.8 

- 

Radius  of  Gyration  for  Pitch 

16.9 

18.3 

Radius  of  Gyration  for  Roll 

10.2 

5.3 

Waterplane  Area,  m2 

193.9 

- 

Maximum  Sectional  Area,  Demihuli,  m2 

- 

40.8 

Projected  Side  Area,  Demihuli,  m2 

- 

- 

Neutral  Axis  Height  from  Waterline,  m 

- 

- 

Heave  Period,  sec 

9.2 

- 

Roll  Period,  sec 

18.7 

6.35 

Model  Scale  Ratio 

1/22.5 

1/11.682 

Once  and  are  known,  <J>B  can  be  obtained  by  Equation  (9) 
and  H  by  Equation  (4).  The  motion  amplitudes  | C j I  which  are  nor¬ 
malized  by  appropriate  factors  are  shown  in  Figure  2  for  beam 
(y  =  90°) ,  bow-quartering  (y  =  135°),  and  head  (y  =  180°)  waves. 

The  ship  length  L  used  in  the  figures  and  text  is  the  length  of 
the  submerged  main  hull,  and  the  notation  S  used  in  the  normal¬ 
ization  of  the  roll  amplitudes  means  one-half  of  the  separation 
distance  between  the  longitudinal  centerplanes  of  individual  hulls. 
It  can  be  noted  in  Figure  2  that  the  non-dimensionalized  peak  heave 
amplitude  occurs  at  the  wave  length  (X)  which  is  about  twice  the 
ship  length.  The  corresponding  peak  values  for  monohulls  of  equiva 
lent  displacement  would  occur  at  about  X/L  =  1.0  to  1.5.  The 
reason  for  the  SWATH  having  its  peak  value  at  longer  wave  lengths 
compared  to  monohulls  is  due  to  its  longer  natural  heave  period 
owing  to  its  small  waterplane  area.  The  large  sway  motion  in 
beam  waves  can  be  observed  at  about  wave  length  being  equal  to 
0.4  times  the  ship  length  which  corresponds  to  about  1.4  times 


the  separation  distance  between  the  twin  struts.  From  the  inter¬ 
mediate  results,  it  was  found  that  in  the  vicinity  of  X  =  0.4L  the 
sway  added  mass  becomes  negative.  This  is  a  particular  feature 
associated  with  twin-hull  ships  indicating  a  ctrong  hydrodynamic 
interaction  between  the  two  hulls.  The  negative  sway  added  mass 
ras  confirmed  for  a  semi-submersible  twin  cylinder  by  experiment 
Ill] .  It  can  be  seen  in  Figure  3  that  the  spiked  beh;  of  the 

sway  motion  at  X  =  0.4L  influences  the  behavior  of  the  steady  mean 
sway  force  and  roll  moment  in  a  similar  manner. 


0.0  0.4  0,8  1.2  1.6  2.0  2.4 


Wave  Lengtb/Sbip  Length  (A/L) 

Figure  2  -  Motion  Amplitude  for  Beam,  Bow-Quartering  and  Head  Waves 
for  SWATH 


Steady  mean  forces  and  moments  on  the  SWATH  subject  to  beam, 
bow-quartering  and  head  waves  are  shown  in  Figure  3.  The  forces 
and  moments  are  normalized  by  the  factors  pgc^L  and  pgc*L2 ,  res¬ 
pectively  . 

As  expected  from  Equation  (5),  the  steady  mean  surge  force 
induced  by  beam  waves  is  negligibly  small,  since  the  second  term 
in  the  right-hand  side  of  the  equation  which  yields  a  dominant 
contribution  is  zero  for  u  *  90°.  The  first  term  involving  the 
square  of  the  Kochin  function  is  usually  of  higher  order  than  the 
second  term  because  H  *  O(^)  for  a  slender  body.  One  can  see  in 
Figure  3  that  the  magnitude  of  the  first  term  of  Equation  (5) 
contributing  to  the  steady  mean  surge  force  is  represented  by  the 
solid  curve  (y  =  90°)  and  is  negligibly  small.  However,  when  the 
wave  heading  is  increased  to  135°  and  to  180°,  the  magnitude  of  the 
mean  surge  force  is  increased.  The  direction  of  the  steady  mean 
surge  force  is  in  the  direction  to  drift  the  ship  backward.  As  ex¬ 
pected,  the  steady  mean  sway  force,  which  will  be  referred  to  as 
drift  force  hereafter,  decreases  as  the  wave  heading  is  changed 
from  the  beam  waves  to  the  bow-quartering  waves.  At  u  *  180°,  the 
drift  force  should  be  negligible  due  to  the  symmetry  of  the  body. 

As  already  explained  earlier,  the  contributions  of  the  term 
involving  | H | 2  to  the  steady  mean  forces  and  moments  are  relatively 
small  compared  to  the  term  involving  the  first  power  of  H.  It  was 
found  from  the  computation  of  ¥l  and  F2  that  the  square  term  of  H 
is  less  than  ten  percent  of  the  first-power  term  of  H.  Based  on 
this  finding,  the  computation  for  F^  is  carried  out  by  neglecting 
f/BB)  since  y^(BB)  as  given  by  Equation  (8b)  is  of  0(Ha).  It  was 
felt  that  the  computation  of  (BB) ^  would  require  two  to 

three  times  greater  computer  time  than  that  required  to  evaluate 

_  /  TO  \ 

F*  ,  does  not  justify  its  merit.  Thus,  the  results  shown  in 
Figure  3  for  the  steady  mean  roll  moment  are  represented  by  F,/10) 
only.  The  change  in  sign  of  the  steady  roll  moment  for  shorter 
wave  lengths  is  worth  noting.  The  positive  roll  moment  is  the 
moment  to  heel  the  ship  toward  the  lee  side  of  the  wave  incidence. 

The  steady  mean  yaw  moment  is  also  shown  in  Figure  3.  For 
beam  waves,  the  steady  mean  yaw  moment  seems  insignificant.  The 
mean  yaw  moment  for  the  bow-quartering  waves  is  found  to  be  the 
same  order  of  magnitude  as  that  of  a  monohull  ship  of  equivalent 
displacement  (2] . 


At  about  X/L  *  0.3  the  steady  mean  forces  and  moments  shown 
in  Figure  3  exhibit  a  singular  behavior.  The  wave  length  corres¬ 
ponding  to  about  0.3L  happens  to  be  equal  to  the  separation  dis¬ 
tance  between  the  inner  sides  of  the  vertical  struts  of  the  SWATH. 

It  is  mathematically  well-known  that  there  exists  a  breakdown  of 
the  velocity-potential  solution  at  the  wave  lengths  which  are 
integer  multiples  of  the  separation  between  two  vertical  barriers. 
This  phenomenon  was  numerically  and  experimentally  confirmed  by 
Lee,  et  al.  [9]  for  infinitely  long  horizontal  twin  cylinders.  In 
the  present  computation,  the  strip  approximation  based  on  the  two- 
dimensional  solution  of  the  velocity  potential  function  has  been 
utilized.  Thus,  the  large  spiked  behaviors  of  the  steady  forces 
at  this  wave  length  are  regarded  to  be  a  direct  consequence  of 
the  aforementioned  shortcoming  of  the  strip  approximation.  Al¬ 
though  the  phenomenon  could  have  been  exaggerated  by  the  two- 
dimensional  approximation,  the  existence  of  such  a  phenomenon  for 
a  SWATH  having  long  surface-piercing  struts  cannot  be  ruled  out 
entirely.  At  DTNSRDC  peculiar  phenomena  have  often  been  observed 
in  the  roll  behavior  of  SWATH  models  in  waves.  These  phenomena 
are  the  asymmetric  roll  motion  and  the  resonant- like  motion 
at  the  wave  lengths  far  shorter  than  the  expected  resonant  wave 
length.  These  phenomena  could  be  related  to  the  rapidly  oscillat¬ 
ing  behaviors  shown  in  Figure  3  for  X/L  <  0.6. 

To  investigate  the  cause  for  the  oscillatory  behaviors  of  the 
curves  in  Figure  3  for  X/L  <  0.6,  computations  are  made  by  setting 
the  motion  of  the  body  to  be  zero,  i.e.,  |?j|  =0  for  j  *  1,  2, 

.  .  .  ,6.  In  this  way,  the  effect  of  the  motion  on  the  steady 
mean  forces  and  moments  could  be  investigated.  Figure  4  shows  the 
steady  mean  sway  force  and  roll  moment  for  restrained  and  unre¬ 
strained  conditions  for  beam  waves.  One  can  observe  that  the  large 
hump  and  hollow  for  the  unrestrained  condition  diminishes  for  the 
restrained  condition  for  X/L  <  0.6.  This  implies  that  for  X/L  <  0.6 
the  motion-associated  hydrodynamic  effects  have  greater  influence 
on  the  steady  forces  and  moments  than  the  wave  diffraction  effect 
does.  That  is,  the  disturbance  potentials  for  j  *  1,  2,  .  .  .  ,6 
have  greater  influence  than  $7  in  shorter  wave  lengths.  This  impli¬ 
cation  can  be  related  to  the  hump-and-hollow  behavior  of  sway  and 
roll  motion  for  X/L  <  0.6  shown  in  Figure  2.  It  is  plausible  that 
the  hydrodynamic  interactions  between  the  two  hulls  which  are  caused 
by  the  sway  and  roll  motion  could  have  a  direct  effect  on  the  oscil¬ 
latory  behavior  in  the  steady  forces  and  moments  in  shorter  wave 
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lengths.  However,  it  can  be  realized  readily  that  for  0.6  <  X/L  < 
1.6  the  trends  and  the  magnitudes  of  the  two  curves  in  Figure  4  are 
not  so  different.  This  implies  that  the  wave  diffraction  effect 
which  entirely  governs  the  steady  mean  forces  and  moments  on  a  re¬ 
strained  body  seems  the  dominant  factor  in  this  wave-length  range. 
As  wave  length  increases,  the  diffraction  effect  should  diminish, 
and  it  seems  from  Figure  4  that  such  a  trend  holds  true. 


W*v»  Length/ Ship  Length  U/U 


Figure  4  -  Steady  Mean  Sway  Force  and  Roll  Moment  for  Restrained 
and  Unrestrained  Condition  for  SWATH 


To  assess  the  maximum  magnitudes  of  the  steady  mean  forces  and 
moments,  a  wave  length-to-amplitude  ratio  of  60  is  assumed,  and  the 
calculated  values  are  shown  in  Table  2  in  which  the  forces  and 
moments  are  respectively  normalized  by  the  displacement  of  the  ship 
and  the  displacement  times  one-quarter  length  of  the  ship.  The 
values  in  Table  2  are  obtained  by  taking  the  values  from  Figure  3. 
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Table  2 

Maximum  Steady  Mean  Forces  and  Momenta  for  A/(a  =  60  for  SWATH 


IF,  I/A 

|F2|/A 

IM/<AV 

|Fel/<A*» 

0.75% 

3.38% 

2.92% 

2.17% 

A/L 

1.9 

1.4 

1.4 

1.6 

P 

135° 

90° 

90° 

135° 

The  heeling  angle  induced  by  the  steady  mean  roll  moment  given 
in  Table  2  can  be  obtained  by 
180 


0  25  -  X 

AGM 


10.6 


where  A  =  2802  long  tons  and  GM  *  2.9  m  (see  Table  1).  The  heeling 
angle  is  toward  the  lee  side  of  the  wave  incidence.  Such  a  steady 
heeling  moment  induced  by  the  body  and  wave  motions  may  explain  the 
observed  asymmetric  roll  motion  of  SWATH  models  in  waves.  That 
is,  when  a  SWATH  model  is  subject  to  regular  beam  waves  at  zero 
speed,  the  oscillatory  roll  motion  often  takes  place  about  the 
vertical  plane  which  is  somewhat  inclined  toward  the  lee  side  of 
the  wave  incidence. 


In  Figure  5  the  steady  mean  sway  forces  of  SWATH  and  a  mono¬ 
hull  having  the  particular  dimensions  as  shown  in  Table  1  are  com¬ 
pared.  Both  ships  have  the  same  length.  To  compare  the  drift 
force  in  a  meaningful  way,  the  drift  forces  are  normalized  by  the 
factor  A(i;a/L)2  and  are  shown  versus  the  ratio  of  wave  length  to 
draft  (T)  of  individual  ships.  It  can  be  observed  that  at  about 
A/T  =  6  two  curves  intersect  and  for  the  smaller  values  of  A/T  the 
monohull  curve  shows  greater  values  than  the  SWATH  curve,  and  the 
trend  reverses  for  the  greater  values  of  A/T.  Since  the  drift  force 
is  proportional  to  the  square  of  the  wave  height,  and  the  wave 
heights  for  longer  waves  would  be  larger  than  for  shorter  wave 
lengths,  it  is  anticipated  that  SWATH  ships  can  be  subject  to 
larger  drift  forces  than  monohull  ships.  This  would  mean  that  if 
a  SWATH  type  of  configuration  is  considered  for  a  drilling  ship, 
a  more  powerful  position  control  device  than  for  an  equivalent 
monohull  ship  may  be  necessary. 
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Figure  5  -  Steady  Mean  Sway  Force  of  SWATH  and  Monohull  Ship 
in  Beam  Waves 


4.  Conclusion 

Within  the  limit  of  the  present  investigation  it  is  found  that 
SWATH  ships  can  be  subject  to  a  larger  drift  force  than  equivalent- 
displacement  monohull  ships  for  wave  lengths  greater  than  six  times 
the  draft  of  the  SWATH.  This  means  that  a  SWATH  may  require  more 
power  for  position  control  than  a  monohull  of  equivalent  displace¬ 
ment.  The  magnitude  of  the  steady  heel  moment  exerted  on  a  SWATH 
ship  could  be  sufficiently  large  to  cause  the  ship  to  assume  a 
heeled  mean  position  resulting  in  an  asymmetric  roll  oscilla¬ 
tion.  The  results  obtained  in  the  present  investigation  indicate 
that  the  wave-diffraction  has  a  significant  effect  on  the  steady 
mean  forces  and  moments  on  a  SWATH  ship  for  2.0  <  A/(2S)  <  5.0. 
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Appendix  I  -  Derivation  of  Mean  Roll  Moment 
The  rate  of  change  of  angular  momentum  of  the  fluid  in  the 
volume  V  enclosed  by  the  surface  S  is  given  by 


x  V  dv  »  p  fff  T 

3V 

*  It  dv 

+  off  L  *  Y  Un  ds 

V 

V 

8 

-»///- 

v 

j_(V  •  V)V  -  ^  -  V(gz)}  dv  + 

»// I  <  v 

S 

•  -’///  [’  • 
y 

*  )V(r  x  V)}  +  r  x 

V  (-J-  +  gz)j  dv 

+  p  //  i  > 
s 

(  V  Un  ds 

■  -p//[<I  > 

■  ?><v»  -  V  *  i 

(r  x  n)  + 

gzr  x  ( 0 , 0 , n , )J  ds  (I 

where  V  is  the  velocity  vector  of  the  fluid,  U  is  the  normal  velo 

n 

city  of  the  surface  S,  VQ  the  normal  fluid  velocity  on  S,  r  the 
position  vector  and  n,  the  z-component  of  the  unit  normal  vector  n 
on  S. 


If  we  let  S  *  S6  +  SF  ♦  where  SQ  represents  the  body  sur¬ 
face,  SF  the  2—0  plane  excluding  the  waterplane  surface  of  the 
body  and  the  side  and  the  bottom  surface  of  a  vertical  cylin 
der  having  a  largo  radius  R,  we  can  show  from  Equation  ( I — 1 )  that 
the  moment  vector  M  about  the  origin  is  given  by 


p  r  x  n  ds 


*  - «//  z(yn, ,-xn, ,0)ds 

S. 

P  (x.y.O)  *  dxdy  -  JJ  (x.y.O)  *  n  p  dx  dy 

Sp 

P //  P<I  *  fi)ds  -  off  (r  *  7*)*n  ds 


(1-2) 


where  V*  »  V, 


and  SR  is  the  side  surface  of  the  vertical  cylinder. 


1) 
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The  x-component  of  the  time-average  of  Equation  (1-2)  is 
obtained  as 


=  -P 


(K  -  t  I  v*i  *)y  dxdy 

-pff  {f  iy*i  2n2  +  (y*z  -  zV*n} ds 


(1-3) 


>n  principle,^/*  yzns  ds  *  0 
>2),  and^J?t#y  dx  dy  * 


since  dK/dt  is  zero  due  to  the  conservation 
due  to  the  symmetry  of  the  body  within  0($a 
{flt  zds  *  0  +  0( $2 ) . 

Using  the  cylindrical  coordinate  system,  x  =  R cos  0  and 
y  =  R  sin  6,  we  can  rewrite  Equation  (1-3)  as 

f„  =  -5-  JJ  y(**  +  **  -  *£>  dxdy 
SF 

-prJ^  d  bJ  n  2 1 7#  I  2  +  (y»2  -  z*y>*n}dz 

Since  the  asymptotic  behavior  of  •  as  R-*°°  can  be  expressed  as 

K7 

♦  -  e  4(x,y,0,t),  we  can  show  that 


F*  =  i//y  <*x  +  *y  '  *z>  dx 


dy 


i,f  n,  (*2  +  *2  ♦  K2*2)  dt 
8K2  Jn  x  y 

CR 


-of  (i  ‘  *  5*  *y)  *■  " 

CR 

where  f  dt  ■  R  de  and  the  elementary  integrals  such  as 

jf “  e2kzdz  =  fa  and  f_l  ze2kz  dz  =  ^  were  used. 


(1-4) 


From  the  Gauss  theorem,  we  have 


L n 


2  A  A  2  x  If  2  a  2  \AQ.  =  *// 


,(*2  ♦  •;  +  K242 )dt 


(«  *  +  #  ♦  ♦  K2**,, )dx  dy 

'  x  xy  y  yy  y 


f  n2(*2  >  4*  +  K2*2)dl, 


(1-5) 
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and 


L  (i*+  5F*y)*n  dt  m  T  ff  {(**y  +  yl9*l*  +  y*9^) 

R  Sp 

*  sh  (’*»  ’  ’♦  *  V '*)}'*=  *Xw  (*♦  *  jp»y)  ♦„  «  «-«> 

where  j£  di  is  the  integral  along  the  contours  of  the  waterplane 
area  of  the  body  in  the  counterclockwise  direction,  and  $  * 

(4-  &)• 

Substitution  of  Equations  (1-5)  and  (1-6)  into  Equation  (1-4) 
yields 

T.  •  -  t//  (»x*  *  V  *  «■»)  (iP  %  *  »♦)  <”"iy 

SF 

*  sih/  *•(' 


*  -*■  *2  K242)d* 


♦  2  +  + 
x  y 


-p/r  (**  *  JP  *»)  *»'“ 


(1-7) 


We  arbitrarily  assume  that  the  velocity  potential  *  can  be 
analytically  extended  onto  the  waterplane  area  S^.  Then,  we  can 
write 


».---§•//  (•***  V  *  K!’)(5P  %  *  ,v) 


iff  K* 


f  n  (♦*  t  ♦’  +  K2*2) 

8K2  Jq  *' *  y  ' 


dl 


d  * +  4K2  *y)  *n 


dt 


(1-8) 


If  we  apply  the  same  procedure  of  deriving  Equation  (1-7) 
from  (1-4)  on  the  domain  Sw,  we  cr.n  show  that 
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£//(*xx  +  % y  +  K’*)(dh  *y  +  y*)  dXdy 

sw 

-  t  //  y  K  +  *y  ■  **)  dx  dy 

sw 

-  £  f  +  k2*2)  de 

-of  (  “o’  ♦  +  — — -  ♦  \  ♦  dl 
PJr'2  4K2  y/  n 


Substituting  the  foregoing  identity  into  Equation  (1-8),  we 
can  derive 

T-  =  -  T  //  (*xx  +  *yy  +  ^(i?  *y  +  *)  dxdy 
SF+SW 

+  T  //  y  (*x  +  *y  '  *z)  dxdy 

sw 

= - ff  <♦*  +  2K2y$*)(92  +  K2)*  dxdy 

8k2  jj  y 

SF+SW 

+  \  ff  y  +  *y*y  -  ^z^z^x  dy  (1-9) 


We  impose  the  atmospheric  pressure  condition  for  on  S^, 
i.e.,  (J>B  =  0  on  and  obtain 


Ft, 


—  ff  <♦*  +  2K2y$*)(92  +  K2)*  dxdy 

8K2  JJ  y 
Vsw 


-  T  //  yl<*Iz  +  *Bz)(*Iz  +  *Bz} 

sw 

-  +  ♦ly^y>idxdy’ 

since  tB  *  0  on  Sw  implies  that  $Bx  «  ♦gy  *  0  on  Sw- 
We  can  show  that 


(1-10) 
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II  y(*Ix*lJ  +  *Iy*lJ  -  ♦lz*li)dxdy 

SW 

*  ff  -  K*>(^)1  dxdy 

SW 

*  -2gKCAa  JJ  y  dxdy  -0 


due  to  the  symmetry  of  the  waterplane  area  of  the  twin  hulls. 
Thus,  Equation  (1-10)  reduces  to 


F*  “  ‘  Sk 1  ff  (*y  *  2K*y**)(9*  *  K*>*  dxdy 

SF+SW 

-*//  yl*Bzl  *  dxdy  -  T  Reff  y  *Iz*Bz  dx  dy 


(1-11) 


If  we  assume  that  4gz  and  $Iz  are  approximately  only  functions 
of  x  on  each  waterplane  of  the  two  struts,  we  can  neglect  the  last 
two  integrals  in  Equation  (1-11).  Furthermore,  since  ($2+K*)$j  ■  0, 
we  finally  obtain 


f  =  -  -£_ 
"  8K2 


ff  ^Uiy  ♦  V* +  2K*y(*i  +  vf^2  +  r2)*b  ^oy  ci-12) 

Sw+S*. 


7  "I 


The  foregoing  expression  can  be  divided  into  two  parts  accord¬ 
ing  to  the  products  of  ^  and  $B  in  the  integrand.  That  is, 

F-(IB>  m  '  8K iff  (*Iy  +  2Kly*J)(91  +  k*)*b  dxdy  (1-13) 


F%<BB)  *  ‘  8K2  //(* By  +  2K2y^)(92  +  K2)*0  dxdy 


(1-14) 


where  Sp  -  Sp  +  Sw. 

The  expression  for  ♦g  derived  from  Green's  theorem  can  be  given 
by  (see,  e.g.,  Equation  (74)  of  15]) 

♦B<*.y.°)  -  ^ f2"**f~  tott*~ik(xcoa*  *  y sln ® )H(k , 0 )  (i-i5) 

where  the  symbol means  that  the  path  of  the  integral  is  taken 
below  the  pole  k«K. 
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where  6  is  the  Dirac  delta  function  and 

l  *  k  cos  6  ,  m  *  k  sin  6 

l'  ■  k'  cos  0  1 ,  m*  »  k*  sin  0 

Then,  we  can  derive  that 
y»(IB)  =  ^^[-sinw  H(K, it-u ) 

+  (sln  0  h  +  nr5  A)  <K+*)H<k,e)|k=K  ] 

0*Tt-u 


Re//  (ZSinw  +  iy)(*Bn  -  *B  £i) 


eK(z  -  ix  cos  u  +  iysinu)dg 


(1-18) 


?><bb>  ■  ^r'd9s‘"9/o'’dk ^ 

u  si? 

'  ImI*’ “h  dk  w  H'(sln9i*  * 

-J^U 

■  55Hf  fo'M sln 6/0" “ 

-^/„2'd9sl"9/„"dk«  |H|1 

-  tst/02'  d9  «*«■•>//  <**‘"9  *  ‘V)(*B„  -  ♦„  £„) 


k* (k+K) 


|  H(k , 0 ) | 


.,K(z  +  ix  cos  0  ♦  iysin0) 


-  iifi/02"  d9/‘ dk  ^  “*//  <dSl"9  *  -  *.fj 


sk(z  +  ix  cos  0  +  iysin0) 


The  first  two  terms  of  the  foregoing  expression  involved  with 
|H| 2  are  pure  imaginary  values;  hence,  they  do  not  contribute  to 
v  }  which  should  be  a  real  function.  Thus,  we  have 


21 


VBB)  -  -  £  ■*//’<*  !*«.•)//  C*..  *  1V)(»B„  -  *B  &) 

So 

eK(z  +  ix  cos  0  +  iysin0)ds 

'  16s*  "*  lC*£*  TT  ■*//  <*s1"9  *  ly> 

So 

/  3  \  k(z  +  ix  cos  0  +  iy  sin  0 )  .  <'1-19 

VBn  '  *B  Snl  e  ds  (I  19 
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